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ABSTRACT 

The  effects  of  uncertainties  in  the  geophysical  and  geodetic  parameters 
related  to  the  establishment  of  ballistic  missile  launch  sites  on  the  impact 
accuracy  of  the  ballistic  missile  are  derived  analytically.  Curves  of  the 
derived  impact  miss  coefficients  are  presented  for  radio-  and  inertially- 
guided  missiles  for  ranges  of  500  to  21,000  nautical  miles  and  burnout 
angles  of  6  to  90  degrees  from  the  vertical. 
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I.  INTRODUCTION 

This  report  documents  a  study  of  those  geophysical  and  geodetic  (GfcG) 
phenomena  which  are  pertinent  to  the  establishment  of  a  missile  launch  site. 
A  study  of  this  particular  class  of  GRtG  effects  is  necessary  for  ascertaining 
their  relationship  to  the  ultimate  accuracy  of  the  ballistic  weapon  system. 

In  particular,  this  relationship  is  relevant  and  valuable  in  the  study  of  pos¬ 
sible  constraints  on  future  weapon  systems.  This  is  not  to  under  estimate 
the  importance  of  these  effects  for  present  weapon  systems;  however,  the 
determination  of  G&G  launch  area  effects  for  presently  well-defined  systems 
can,  and  has  been,  adequately  handled  by  direct  computer  simulation.  The 
study  of  GfcG  effects  for  future  weapon  systems  cannot  be  so  easily  simulated 
since  much  of  the  relevant  information  required  by  the  computer  is  unspeci¬ 
fied.  The  purpose  of  this  report,  then,  will  be  to  provide  a  study  of  the 
particular  class  of  G&G  effects  associated  with  launch  site  development  in 
sufficient  generality  to  allow  reasonably  accurate  estimates  to  be  obtained 
for  any  missile  system  which  might  be  considered.  What  is  gained  in 
generality  will  be,  of  course,  lost  in  accuracy.  For  many  studies  of  future 
systems,  estimates  of  errors  with  an  accuracy  greater  than  ilO  per  cent  is 
inappropriate  and  unwarranted;  this  guideline  was  maintained  in  conducting 
this  study. 

The  GfcG  information  which  must  be  provided  prior  to  the  liftoff  of  a 
ballistic  missile  is  the  location  of  the  launch  site  and  the  orientation  in  space 
of  some  nominal  coordinate  system.  Consider,  specifically,  a  right-handed 
orthogonal  coordinate  system  whose  origin  is  at  the  intended  launch  point; 
define  one  axis,  the  y  axis,  of  this  coordinate  system  as  parallel  to  the 
radius  vector  from  the  center  of  the  earth  to  the  launch  point;  define  the  x 
axis  as  parallel  to  the  plane  defined  by  the  intended  velocity  vector  at  burnout 
and  the  center  of  the  earth;  the  z  axis  is  thus  defined  by  the  requirement  that 
the  system  be  right-handed.  This  coordinate  system  just  defined  shall  be 
termed  the  nominal  coordinate  system.  Several  unpleasant  properties  of 
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auch  a  nominal  coordinate  system  should  be  mentioned:  (1)  except  at  special 
points  (the  poles  and  equator)  the  y  axis  is  not  normal  to  any  reference 
datum  surfaces  presently  in  use;  and  (2)  most  ballistic  missile  trajectories 
in  practice  are  not  coplanar  due  to  the  earth's  rotation  thus  the  burnout 
velocity  vector  does  not  lie  in  the  xy  plane. 

Properties  1  and  2  above  result  from  comparing  an  oblate,  rotating 
earth  (ORE)  with  a  spherical,  nonrotating  earth  (SNRE) .  A  key  assumption 
of  this  report  will  be  that  deviations  caused  by  G&G  effects  from  a  trajectory 
traversed  over  an  ORE  can  be  approximated  by  similarly  caused  deviations 
from  a  trajectory  traversed  over  a  SNRE.  A  qualitative  argviment  in  sup¬ 
port  of  this  assumption  may  be  made  as  follows;  for  ICBM's  the  differences 
between  trajectories  over  an  ORE  and  a  SNRE  are  small  as  compared  to  the 
trajectories  themselves;  consequently,  the  deviations  caused  by  small  G&G 
effects  for  SNRE  trajectories  will  differ  from  similarly  generated  deviations 
for  ORE  trajectories  by  a  small  perturbation  in  the  already  comparatively 
small  difference  between  ORE  and  SNRE  trajectories.  The  assumption  that 
it  is  sufficient  to  treat  SNRE  trajectories  removes  the  necessity  of  considering 
the  disparities  indicated  by  properties  1  and  2. 

If  the  ICBM  is  considered  to  follow  a  specified  trajectory  in  the  nominal 
coordinate  system,  then  the  class  of  G&G  uncertainties  under  consideration 
have,  at  liftoff,  the  effect  of  rigidly  transforming  the  coordinate  system  in 
which  the  missile  actually  flies.  There  are  six  and  only  six  independent 
transformations  of  the  nominal  coordinate  system:  three  rotations  about 
each  of  the  axes;  three  translations  along  each  of  the  axes.  Each  of  these 
transformations  can  be  identified  as  being  generated  b/  a  particular  G&G 
uncertainty  as  shown  in  Table  1. 
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Transformation 
Rotation  about  z  axis 

Translation  along  y  axis 

Translation  along  x  axis 

Rotation  about  x  axis 

Rotation  about  y  axis 
Translation  along  z  axis 


Table  1. 

Corresponding  G8tG  Uncertainty 


Deflection  of  the  vertical  uncertainty; 
component  in  trajectory  plane 

Elevation  above  reference  ellipsoid 
uncertainty 

Horizontal  position  uncertainty;  down- 
range  component 

Deflection  of  the  vertical  uncertainty; 
component  perpendicular  to  trajectory 
plane 

Azimuth  uncertainty 

Horizontal  position  uncertainty;  cross¬ 
range  component 


A  complete  evaluation  of  the  effects  of  the  above  uncertainties  requires 
that  they  be  interpreted  in  terms  of  target  miss  of  the  missile.  The  required 
analysis  may,  for  simplicity,  be  divided  into  two  phases:  Phase  1,  the  analy¬ 
sis  of  the  G&G  effects  through  the  powered  portion  of  the  missile  flight  to 
burnout  at  which  the  perturbation  of  the  burnout  conditions  should  be  obtained; 
Phase  II,  the  analysis  of  the  effect  of  previously  defined  burnout  perturba¬ 
tions  on  the  free  flight  of  the  missile  until  impact.  Re-entry  conditions  will 
be  ignored  under  the  assvimption  that  deviations  of  the  missile  from  the 
nominal  re-entry  point  will  produce  an  equal  error  in  impact  point.  The 
analysis  in  Phase  I  will  be  considerably  simplified  by  use  of  the  facts  that 
the  geocentric  angle  between  liftoff  and  burnout  for  most  missiles  is  of  the 
order  of  5  degrees  and  the  burnout  height  is  generally  of  the  order  of  1/10, 
or  less,  of  the  earth's  radius  and  consequently  both  quantities  may  frequently 
be  ignored. 
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II.  PHASE  I  ANALYSIS 

From  liftoff  to  burnout  a  missile  is  acted  upon  by  the  thrust  force  of 
the  rocket  engines,  aerodynamic  forces,  and  gravitational  forces.  The 
gravitational  forces  are  wholly  a  function  of  position  of  the  missile.  The 
aerodynamic  forces  are  strongly  a  function  of  time  and  only  weakly  a 
function  of  position.  The  thrust  force  can  be  a  more  complex  function  of 
time  and  position,  since  this  force  is  controlled  by  the  guidance  system. 
Two  basic  categories  of  thrust  guidance  may  be  discriminated;  (1)  the 
thrust  is  controlled  so  as  to  conform  to  a  predetermined  acceleration 
program:  such  a  process  would  be  characteristic  of  an  all-inertial  guid¬ 
ance  system  and  the  thrust  force  is  almost  entirely  a  function  of  time; 

(2)  the  thrust  is  controlled  in  a  manner  which  will  maintain  the  missile  on 
a  predetermined  trajectory;  such  a  process  requires  a  guidance  system 
which  can  provide  feedback  control  to  the  missiles,  a  process  which  is 
characteristic,  ideally,  of  radio  guidance  systems  and  partly  character¬ 
istic  of  astro-inertial  and  radio -inertial  guidance  systems.  Since  it  will 
be  shown  that  the  analysis  for  missile  systems  of  Category  2  may  be 
obtained  by  a  specialization  of  the  results  for  Category  1  missile  systems, 
the  Category  1  missile  systems  will  be  considered  first. 

A.  Category  1 

If  the  position  vector  to  the  missile  in  the  nominal  coordinate 
system  is  denoted  by  r,  and  time  differentiation  by  dot's  over  the  letter, 
then  Category  1  systems  involve  the  consideration  of  equations  of  the 
nature  of  (1). 

r  =  F(t)  +  G(r)  r(0)  =  0,  r(0)  =  0  (1) 

where  F(t)  represents  the  thrust  and  aerodynamic  forces  and  G(r) 
represents  the  gravitational  forces. 
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1 .  Deflection  of  the  Vertical,  Downrange  Component  6 1| 

Errors  are  created  by  this  effect  in  the  following  manner: 
the  angle  a,  between  the  thrust  vector,  whose  magnitude  is  F(t),  and  the 
y  axis  is  a  function  of  time  chosen  so  as  to  provide  the  proper  burnout 
conditions;  however,  the  angle  Cl  is  actually  measured  throughout  the  flight 
with  reference  to  a  coordinate  system  rotated  about  the  axis  by  6  j|  ,  see 
Figure  1. 


Figure  1. 

The  thrust  components  in  the  nominal  system  are  thus: 

F^  =  sin  (a  -  6j|)  F(t),  Fy  =  cos  (a  -  6||)  F(t)  (2) 

but  F°  =  sin  CX  F(t),  F°  =  cos  a  F(t)  where  F°,  F°  are  the  components 
which  were  intended  to  be  obtained,  the  nominal  components.  Equation  (2) 
becomes  then 

F^  =  cos  6|[  F°  -  sin  6||  F°,  F^  =  cos  6||  F°  +  sin  6||  F°  (3) 
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Note  that  use  has  been  made  of  a  prior  assumption  that  the  xy  plane 
contains  the  plane  of  the  trajectory.  For  the  sake  of  brevity,  define  the 
operators  T,  as: 

T  T  t 

TCf(t)]  =  J  f(t)dt.  T^[f(t)]  =  /  dt  J  f(§)  d§  . 

o  o  o 

The  operator  T  operating  on  any  constant,  c,  gives  T  [c]  =  c  T  and 
2  2 

T  [c  ]  =  c  (T  /2),  where  T  is  the  total  time  of  integration;  in  this 
report  T  equals  the  time  from  liftoff  to  burnout  of  a  missile.  Then 
from  Equation  (1),  utilizing  x(0)  =  y(0)  =0,  and  representing  a  vari¬ 
able's  value  at  burnout  by  a  subscript  B: 

Xg  =  T^CF^(t)  +  G^(x.y)],  yg  =  T^CF^Ct)  +  Gy(x,y)]  .  (4) 

Combining  (3)  and  (4)  now  gives: 

Xg  =  cos  6||  T^[F°{t)]  -  sin  6j|  T^CF°(t)]  +T^[G^(x,y)] 

(^) 

Yb  =  cos  6|j  T^CF°(t)]  +  sin  6|j  T^CF°(t)]  +  T^[Gy(x,  y)] 

Noting  that  T^[F°(t)]  =  x°  -  T^Cg  (x°,y°)]  , 

X  JD  X 

T^CF°(t)]  =  y°  -  T^CGy(x°,y°)]  (6) 

where  x°,  y°  are  the  nominal  coordinates,  and  anticipating  small  6|| 
so  that  cos  ®||  =  1*  sio  6||  =  6jj  ,  then  (6)  and  (5)  may  be  combined  to 
give: 

Xg  =  ■  *11  ®|| +  ■'■^Cg^(x,  y)] 

(7) 

Yg  =  Vb  ■  y°)]  +  *11 “  ^11 y)^ 


By  changing  the  operator  notation  the  velocity  components  are  similarly 
derived: 
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3  =  *3  -  T[G^(x°,y°)]  -  6||y°  +  &  |jT[Gy(x°,  y°)  ]  +  T  [G^(x,  y)  ] 


'B 


Now 


Yg  -  .y  )]  +  6||X3  -  6||T[G^(x°,  y°)]  +  T[Gy(x,  y)] 


(8) 


^x<**  y'  -  y°)  =  -55r 


I  n  ■'■  ^y  ~X~~  I  r.  “  ^xG°  +  AyG° 

x=x°  ^  x=x° 

o  o 


y=y 


y=y 


and  similarly  G  (x,  y)  -  G  (x*^,  y°)  =  AxG°  +AyG°  where  Ax  =  x  -  x°, 

A  o  ^  ^  o  Q  yy 

^y  =  y  -  y  as  well  as  Ax  =  x-  x,  Ay=y-y  .  Equations  (7)  and  (8)  are 
now  considerably  simplified  to 


Axr  = 

A  O 

-‘ll  ''b 

AyR  = 

A  O 

H’‘b 

= 

A  .  O 

■‘II  ’'b 

AyB  = 

A  •  O 

‘||»B 

^IIYr  +  S||T  [Gy(x°,y°)]  +  t'^CAxG®^  +  ^yG°y] 
6|jT2CGx(x°.y°)]  4  t2[AxG°^  4  AyG^^] 


(9) 


4  6||T  [Gy(x'’,y‘")]  4  t[AxG^  4  AyG°„] 


o  o. 


,o 

Jt 

yx 


,o 

xy 

,o  - 

^yy' 


(10) 


obtained; 


From  Figure  2  the  derivation  for  G^,  G^,  can  be  readily 
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[^1 

^^B  =  ®Pb  ■"  ^go"  1^1 


Now,  from  Figure  2, 

y°  =  (R  +  h°)  cos  i(t°  -  R,  x°  =  (R  +  h°)  sin  i|f°  (14) 

Equations  (12)  and  (13)  are  essentially  integral  equations  which 
converge  quite  rapidly  for  reasonable  values  of  T;  one  iteration  of  (12)  gives; 


6 II  [(R  +  h°)  cos  (i° 


-  *11  go  ^  -  go""^ 


[(l  +  cos  -  1  j  -  6||  ^  [^j  j 

(15) 

B  ^go"^^  *11  (*  Tt)  'I'®  +  j 


Ayg  =  6m  (R  +h°)  sin  +  ^g^T^  6| 


Th'  terms  under  the  T  operator  in  (15)  are  already  of  the  order  of  terms 
previously  neglected,  hence: 


^^B  =  -*||  Tb  +  g 


=  6||R  sin  H-g 


T^l 

o  "Z") 


If  (16)  is  substituted  into  (13)  the  terms  under  the  T  operator  are  clearly 
negligible,  giving: 

=  -*11  ‘^^B  go^^  =  -*ll  <^B  +  '^'b^  ^  go"^^ 

^^B  =  *pb  =  *irB  <^B 
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The  relation  between  6j|  and  perturbations  in  the  burnout 
conditions  are  given  by  (16)  and  (17)  but  in  a  clumsy  reference  system  for 
the  analysis  in  Phase  II.  More  useful  quantities  would  by  Ah,  A3,  Av 
rather  than  Ax,  Ay,  Ax,  Ay;  the  necessary  relationships  can  be  quickly  de¬ 
rived  using  Figure  2: 

tan  t  =  ^Ji_,  d\|f  =  (R  +  h)  ^  [cos  (r  dx  -  sin  #  dy]  (18a) 

2  2  2 

(R  +  h)  =  X  +  (R  +  h)  ,  dh  =  sin  i)  dx  +  cos  fdy  (18b) 
tan  (3  4  (()  =  ^,  d3  +  dt  =  —  [cos  (3  4  f)  dx  -  sin  (3  4  1))  dy]  (18c) 

y  ^ 

=  x^  4  y^,  dv  =  sin  (3  4  t|t)  dx  4  cos  ( 3  4  (t)  dy  (18d) 


{19d) 


2.  Elevation  Uncertainty,  C 

The  effect  of  this  uncertainty  is  essentially  to  translate  the 
origin  of  the  nominal  coordinate  system  a  distance  e  along  the  y  axis. 
The  equations  for  x  and  y  corresponding  to  (4),  using  the  T  operators. 


are 
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Xr  =  T^[F^(t)  4  G^(x.y)].  Yg  =  [Fy(t)  4  Gy(x.  y)]  4  6  (20) 

Noting  that  F  (t)  =F°(t),  F  (t)  =  F°(t)  the  integral  equations  for  Ax,  Ay 
X  X  y  y 

may  be  derived  in  a  fashion  similar  to  that  used  in  obtaining  (9)  and  (10) 
to  obtain 


‘■‘b  =  (21) 

iyg  .  +  iyG“^]  *  « 


Ax  =  T  [AxG°  4  AyG°  ] 
B  XX  '  xy 

Ay„  =  T  [AxG°  4  AyG°  ] 
^  B  yx  ^  yy 


(22) 


Equations  (11)  may  be  used  to  effect  an  immediate  simplification  in  (21) 
and  (22) 


From  (23)  and  (24) 


and  hence 


Axb  = 

llfl 

(23) 

^Yb  = 

Axb  = 

1^1 

(24) 

^Yb  = 

it  is  clear  that 

Ax  i 

0  Ay  i  e 

(25) 

Ax  i  0 

=  2g^  j 

(26) 

Page  12 


Then,  using  Equations  (18a,  b,  c,  d)  with  (25)  and  (26)  gives 


At 

1 

"  R  4  h 

—  [-  sin  tx 

b 

Ah 

=  ®  cos 

♦b  =  ® 

A8 

c 

=  ■  Ti 

sin  (3g  4  \ 

Av 

=  i  +  *b' 

2g  T 


3.  Horizontal  Position  Uncertainty,  Downrange 


Component,  d 

The  effect  of  this  uncertainty  is  essentially  to  translate  the 
origin  of  the  nominal  coordinate  system  a  distance  d  along  the  x  axis. 
The  equations  for  x  and  y  now  become 

Xfi  =  CFx(t)  +  G^(x,y)]  +  d 

Yb  =  +  Gy(x,y)]  (28) 

CF^(t)  +  G^(x,  y)],  Vg  =  T  [F^(t)  +  Gy(x,  y)]. 

Again  noting  that  F^(t)  =  F°(t),  Fy(t)  =  F°(t)  we  obtain: 

Ax„  =  d  4  T^CAxG®  +  AyG°  ]) 

B  XX  '  xy 

(29) 

Ay„  =  [AxG°  4  AyG°  ] 


Ax„  =  T  ^xG°  4  AyG°  ] 
B  XX  '  xy 


Ay_  =  T  l>xG°  4  AyG°  ] 


Using  Equations  (11),  Equations  (Z9)  become 


Equations  (31)  in  (30)  give; 


(31) 


(32) 


Equations  (18a,  b,  c,  d)  may  now  be  used  to  obtain: 

Ahg  =  d  sin  =  0  (33b) 

cos  (8^  +  tg)  j  (33c) 

=  -  ^  goT  sin  (3g  +  i)g)  (33d) 


d  r  ^o 
-  _  1  +  - 


4. 

Analysis  of  this  effect  is  facilitated  by  introducing  the  arti¬ 
fice  of  letting  the  nominal  trajectory  plane  make  an  angle  0  with  xy 
plane  and  later  allowing  9  to  approach  0.  From  Figure  3, 


Deflection  of  the  Vertical  Uncertainty,  Crossrange 
Component,  6j 


Figure  3. 
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F^(t)  =  F(t)  sin  (0  +  6  ) 

F^(t)  =  F(t)  cos  (0  4  6  )  (34) 

F^(t)  =  F“(t) 

Note  that  F(t)  Bin  0  =  F°(t)  and  F(t)  cos  9  =  F®(t)  thus  Equation  (34) 
gives 

F  (t)  =  F°  cos  6j^  4  F°  sin  6  i  =  F°  +  6  ,  F° 

z  z  y  z  ^  y 

F  (t)  =  F°  cos  6j_  -  F°  sin  6j_  =  F°  -  (>±  F° 

y  y  z  y  X 

F  (t)  =  F° 

X  X 

In  the  operator  notation,  we  have: 

^  +  G^(x,y.  z)],  [F^(t)  4  Gy(x,  y,  z)  ] 

+  G^(x.y,z)]  (36) 

plus 

=  x|  -  T2tG^(x°,y°,z‘>)] 

^^[Fp  =  y°  -  T^[Gy(x°,y“,z°)]  (37) 

T^[F°]  =  4  -  T2[G^(x'’,y'>,t'>)] 

Remembering  that  G  (x,  y.  z)  -  G  {x°,y°,z°)  =  AxG°  4  AyG°  4  AzG° 

X  X  XX  '  xy  x; 

etc.  ,  and  using  (37)  and  (35)  in  (36)  Equations  (38)  are  obtained: 

Axg  =  t2[AxG^  +  AyG°y  4  AzG°^] 

^Yb  =  "^J-^B  ^  [G^(x°,  y°,  z°)]  4  T^C^xG^  4  ^7^°^  +  6z  G°y.] 

=  ^-lYb  -  *-lT^fcy<x°.y°,z°)]  4  T^CAxG®^  4  6yG°y  +  AzG°^] 

(38) 
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A  similar  development  replacing  T  with  T  gives 
Ax  =  T  [AxG°  +  AyG°  +  Az  G°  ] 

Ay  =  -6J.Z3  +  [G^(x°.  y°.  z°)]  +  T  +  ^yGyy  + 

y R  -  CG  (x°.  y°,  z°)  ]  +  t[AxG®  +  ^yG°  +  AzG°  ] 

.D  y  ZX  ^  j  ZZ 


Allow  0  to  approach  zero;  now  z°,  z®,  G^(x®,  y®,  z®)  -  0  and  since 


,0 

xz 

,0 

Vz 

GM 

3  GM  X®  z® 


(x®^  +  (R  +yO)^  +  z®^] 

_ 3GM(R  +  y®)  z° 

x^  +  (R  +  y®)^  4  zo^]^ 


= 


3  GM  z‘ 


zz  f  ,  y  f  7  ?  ?l5/2 

lxO'^  +  {R+y)‘^  +  z°  J  |x  +  (R  +  y)  +  7r\ 

we  also  have,  in  addition  to  (11), 

G®  =  G®  =  G®  =  G®  =0,  G®  =  -  ^  (40) 

xz  yz  zy  zx  zz  R  '  ' 

Using  these  results  and  Equations  (11),  Equations  (38)  and  (39)  simplify  to: 


=  -0^^  [ 

Ax] 

ttJ 

i  0 

1^1 

i  0 

=  ^-LVb  + 

go^-L 

"T 

6  o  T  1 

Vb  +  ®o  -rj 


I 
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and  using  (41)  in  (39) 


Lk  =  0 


Ay  =  0 


Az  =  +  6^g^T  -  g^T  [^] 


Moreover  y  =  (R  +  h)  cos  (r  -  R  =  h  so  that 


So  ^  Q 

^  ~r  ~  ° 


and  hence 


Ass  =  6j.  [y°  +  g„T] 


=  6j.  Ihp  +  ^  T^] 


For  small  x,  Az  is  the  crossrange  displacement  at  burnout; 
however,  Az  can  most  readily  be  handled  when  translated  into  azimuth  un¬ 
certainty.  Assume  first  that  the  trajectory  plane  makes  an  angle  A  with 

2S 

the  xy  plane,  and  that  the  velocity  vector  lies  wholly  in  the  trajectory 
plane,  then 


tan  A  =  — 
z  5c 

^r\ 


using 


V  sin  (3  +  i;)  cos  A 


z^  =  V  sin  (3  +  ♦)  sin  A^ 


V  sin  (p  + 


(cos  Adz  -  sin  A  dx  ) 
z  o  z  o 


■ 


(44) 


and  now  letting  =  0 


AA 

z 


Az 

V  sin  (P  +  t) 


The  other  burnout  parameters  remain  essentially  unaltered  by  Az  when 
A^  =  0.  The  AA^  in  (44)  is  the  variation  in  azimuth  at  the  liftoff  point;  for 
the  later  analysis  the  azimuth  perturbation  at  burnout  is  desired.  Con- 
sider  Figure  4:  v  is  the  nominal  velocity  vector  which  lies  wholly  in  the 
xy  plane;  v  is  the  velocity  vector  with  3,  and  magnitude  equal  to  v, 
which  is  however  rotated  from  v  through  the  angle  AA^  about  the  y  axis. 
From  spherical  trigonometry 


and 


sin  AA  sin  AA.,_ 

_ z  _  _ zg 

sin  b  “  sin  TFT^r 

cos  b  =  cos  (3  4  t)  cos  (f  4  sin  (P  4  (i)  sin  (r  cos  AA 

Z 


» 


Figure  4. 
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For  small  AA  ,  cos  b  =  cos  0  and  therefore 


or 


sin  AA  _  =  sin  AA 

z,  B  sin  P  z 


LA  ^  LA 

z,  B  sin  P  2 


(45) 


AA^  g  is  the  desired  perturbation  of  the  azimuth  at  burnout.  Putting  (44) 
in  (45)  gives,  including  the  subscript  B  now  for  conformity  with  the  pre¬ 
vious  results: 


AA 


z,  B 


Now  Yg 


Vpj  cos  (0^  +  i)q),  so  we  have 

_ ^ 

sin  0 


AA  =6 
z 


cos  (0g  +  t) 


go  T 

^B  ^B 


(46) 


(47) 


5.  Azimuth  Uncertainty,  a 

This  uncertainty  manifests  itself  by  rotating  the  trajectory 
plane  from  the  nominal  xy  plane  by  an  angle  a  about  the  y  axis.  The 
thrust  functions  may  be  immediately  written  as 

F^(t)  =  F°(t)  -  aF°(t) 

Fy(t)  =  F°(t) 

Fz(t)  =  F°(t)  +  aF°(t) 


and  by  a  development  similar  to  previous  derivations  in  this  report 

Ax  =  aT^CG°(x,  y,  z)]  +  T^[AxG°  +  AyG°  +  Az  G°  ]  -  az° 
B  z  '  xx'xy  xz  B 

Ay-  =  T^[AxG°  +  AyG°  +  Az  G°  ] 

B  yx  ^  yy  yz 

Az  =  ax°  -  aT2[G°(x,y,z)3  +  T^CaxG®  +  AyG°  +  Az  G°  ] 
D  a  X  '  zx  ^  zv  zz 


(48) 


o 


I  c 
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=  -O-Zo  +  IT  Cg°(x,  y,  z)  ]  +  T^[AxG°  +  AyG°  +  Az  G°  ] 
B  B  z  ■'  XX  ^  xy  xz 

Ayg  =  t[AxG^  +  AyG°y  +  G°  J  (• 

Azg  =  ax°  -  aT[G°<x.y.z)]  +  T  [AxG^^  +  AyG°y  +  Az  G°^] 
Now  using  Equations  (11)  and  (40)  to  simplify  (48)  and  (49) 


=  -g 
®o 

(^1 

=  0 

"Vb 

1^1 

=  0 

Az  1 

rt  J 

=  a(R 

Axb  = 

-g 

®o 

1^1 

i  0 

AVg  = 

2g  T 

1^1 

i  0 

Azg  = 

axg  - 

g 

&0 

1^1 

From  (46)  and  (51) 


AA  „  =  a 
z,  B 


sin  (3  +  'I') 


6.  Horizontal  Position  Uncertainty,  Crossrange 
Component,  c 

This  uncertainty  has  the  effect  of  translating  the  actual  co¬ 
ordinate  system  a  distance  c  from  the  nominal  coordinate  system  along 
the  z  axis.  The  basic  equations  for  the  coordinates  are  thus 

Xg  =  [F°(t)  +  G^(x,  y,  z)] 

Yg  =  [F°(t)  +  Gy(x,  y,  z)] 

z  =  [F°(t)  +  G,(x,  y,  z)]  +  c 
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B.  Category  2 

Category  2  missile  systems,  those  with  feedback  to  maintain  them¬ 
selves  upon  a  given  trajectory,  behave  only  as  if  the  nominal  coordinate 
system  were  rigidly  transformed  by  the  G8iG  effect;  o^her  dynamic  con¬ 
siderations  which  arise  from  the  fact  that  the  missile  experiences  a  dif¬ 
ferent  from  nominal  gravity  field  are  removed  by  the  feedback  in  the 
guidance  systems,  i.  e.  ,  the  missile  is  constrained  to  move  exactly  eilong 
the  nominal  trajectory  in  the  transformed  coordinate  system.  The  results 
of  the  analysis  for  Category  1  missile  systems  can  thus  be  used  for  Cate¬ 
gory  2  systems  by  removing  the  gravity  effects;  letting  g^  =  0. 

The  results  of  the  Phase  I  analysis  are  summarized  in  Table  2. 
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Table  2. 


Geodetic  and 
Geophysical 
Uncertainty 

Burnout 

Pe  rturbation/ 
GfcG  Uncertainty 

Category  1 

Category  2 

II 

0 

0 

Ah 

^11  ■ 

R  sin  t 

R  sin  t 

Deflection  of 

A@ 

TTT  “ 

g  T 

-  [  1  +  -2—  cos  (|B+i(i)] 

-1 

the  Vertical 

Downrange 
Comp.  ,  6 II 

II 

-g^T  sin  (e+\(i) 

0 

AA° 

0 

0 

II 

0 

0 

Elevation 
Uncertainty,  e 


II 

0 

II 

1 

AS 

2  g 

®o 

€ 

vTT 

Av 

T-  = 

2g  T 
®o 

aa° 

c  “ 

0 

Az 

0 

T 

—  sin  (3+ilf) 
cos  (e+*) 


0 

1 

0 

0 

0 

0 


(  ) 


o 
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Table 

2.  (Continued) 

Geodetic  and 

Burnout 

Geophysical 

Uncertainty 

Pe  rturbation/ 
G&G  Uncertainty 

Category  1 

Category  2 

II 

1 

R 

1 

R 

II 

0 

0 

Horizontal 

Position 

II 

g  X 

-  [l  4  cos  (P4f) 

1 

■ 

Downrange 
Comp.  ,  d 

Av 

■a"  = 

g  T 

- sin  (P4i) 

0 

0 

0 

Az 

■a"  = 

0 

0 

At 

0 

0 

Ah 

0 

0 

Deflection  of 
the  Vertical 

AP 

0 

0 

Crossrange 

Component, 

Av 

31  = 

0 

0 

AA° 

cos  { 04ilr)  1  ®o  ^ 

cos  (34'I») 

TT 

sin  g  V  sin  P 

sin  P 

Az 

31  = 

,  ®o  2 

h  4  -j-  T 

h 
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Table  2.  (Continued) 


Geodetic  and 

Burnout 

Geophysical 

Pe  rturbation/ 

Category  1 

Category  2 

Uncertainty 

GStG  Uncertainty 

M _ 

0 

0 

a 

0 

0 

a 

Azimuth 

0 

0 

Uncertainty,  a 

It 

0 

0 

AA° 

sin  (3+il;) 

sin  (34i)i) 

a 

sin  ? 

sin  3 

It 

R  sin  t 

R  sin  f 

0 

0 

C 

Ah 

0 

0 

c 

Horizontal 

0 

0 

Position 

Crossrange 

Component,  c 

"c” 

0 

0 

AA° 

-g  T 

0 

c 

R  V  sin  P 

Lz 

1 

1 

c 

O 
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III.  PHASE  11  ANALYSIS 

A  complete  analysis  of  GfcG  launch  area  effects  requires  that  these 
effects  be  interpretable  in  terms  of  target  miss.  The  Phase  I  analysis 
transformed  the  launch  area  uncertainties  to  uncertainties  in  burnout  con¬ 
ditions.  The  analysis  of  this  section  will  concern  itself  with  transforming 
burnout  uncertainties  to  target  impact  uncertainties.  For  the  Phase  II 
analysis  as  in  the  Phase  I  analysis  the  assumption  is  maintained  that  a 
spherical  nonrotating  earth  is  a  sufficiently  close  approximation  to  reality 
to  provide  the  desired  results  within  the  approximate  10  per  cent  guide¬ 
line  mentioned  earlier. 

The  Phase  I  analysis  found  that  launch  area  G&G  uncertainties  were 
interpretable  in  terms  of  six  quantities;  Ai|/,  the  variation  in  the  geocentric 
downrange  burnout  angle;  Ah,  the  variation  in  the  burnout  height  above  the 
sphere;  Ag,  the  variation  in  the  angle  the  burnout  velocity  vector  makes 
with  the  radius  vector  at  burnout;  Av,  the  variation  in  the  burnout  velocity; 

AA°,  the  variation  in  the  angle  the  plane  of  the  trajectory  makes  with  the 
meridional  plane  at  burnout;  Az,  the  variation  in  the  position  of  the  burnout 
point  measured  perpendicularly  to  the  plane  of  the  trajectory.  The  first 
four  quantities  Ai|r,  Ah,  AP,  Av  all  are  perturbations  in  the  plane  of  the 
nominal  trajectory  and  consequently  for  a  nonrotating  earth  these  pertur¬ 
bations  should  only  produce  impact  errors  in  the  nominal  trajectory  plane. 

The  component  of  impact  error  in  the  nominal  trajectory  plane  is  generally, 
and  will  be  here,  termed  downrange  miss.  The  latter  two  burnout  variations, 
AA°,  Az,  are  burnout  perturbations  out  of  the  trajectory  plane  and  consequently 
can  be  expected  to  produce  impact  errors  out  of  the  nominal  trajectory  plane. 
The  component  of  impact  error  measured  from  the  nominal  trajectory  plane 
is  generally  termed  the  crossrange  error.  The  Phase  II  analysis  thus 
naturally  subdivides  into  the  two  categories:  downrange  miss  and  cross¬ 
range  miss. 
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A.  Downrange  Miss 

The  simplest  approach  to  analyzing  the  inplane  burnout  variations 
turns  out  to  be  the  use  of  the  ballistic  missile  "hit  equation" 


r 

o 


1  -  cos  i> 

- 2 - 

'o'"  .  2  p 

-GW  ^ 


+  (3  - 

sin  B 


where  the  quantities  are  as  defined  in  Figure  5. 


Figure  5. 


(58) 


For  most  ballistic  missiles,  if  h  is  the  burnout  height,  r  =  R  +  h,  and 

o 

r^  =  R.  Let  GM  -  M  for  the  sake  of  convenience  then  (58)  becomes 


(R  +h)  sin  8  1^1  + 


sin  8  -  sin  (8-^)1  =  1  -  cos  ^ 


(59) 
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Since  the  downrange  miss  is  just  miss  coefficients  for  h,  3,  v  may 

immediately  be  obtained  by  implicitly  differentiating  (59)  with  respect  to 
h,  3,  V,  and  •p  and  then  substituting 


2 

(R+h)  ^ 


(60) 


into  the  resulting  expressions  to  remove  the  velocity  dependence.  This 
process  gives; 


"  ^  ®  ‘  col  6 .  i  sin 

R  1  -  cos  p 


“D,  C  R  II  =  R 


->  Q  /i  .  I'  \  sm  (23  -  P) 
2  cos  3(1+-=  )  -  } — a — — 

CO. 

MD  ^  ^ 

V  dv  ”  V  Q  K  ^  sin  p 

cos  P  +  -K  sin  3  -! - rr^-T 

R  1  -  cos  p 


From  Figure  5  it  is  clear  that 


=  R  If  =  R 


(61) 


(62) 


(63) 


(64) 


B.  Crossrange  Miss 

For  a  spherical  nonrotating  earth  the  trace  of  the  missile's  tra¬ 
jectory  on  the  surface  of  the  earth  is  a  great  circle  and  thus  the  problem 
of  determining  crossrange  miss  becomes  one  of  spherical  trigonometry. 
With  no  loss  in  generality  a  spherical  coordinate  system  may  be  adopted 
whose  equatorial  plane  includes  the  nominal  trajectory  plane.  Figure  6 
should  serve  to  clarify  the  particular  coordinate  system  used.  In  Figure  6 
the  point  B  is  the  nominal  burnout  point  projection;  B'  is  the  actual  burnout 
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point  projected  onto  the  sphere;  ^  is  the  geocentric  range  angle  as  in 
Figure  5;  A'  is  the  actual  burnout  azimuth  at  B'.  The  point  I  is  the  nominal 
impact  point  and  I'  is  the  actual  impact  point.  From  the  previous  defini¬ 
tions  A0  is  the  crossrange  miss,  and  AX  is  the  downrange  miss. 


The  cosine  law  of  spherical  trigonometry  gives: 


cos 


(a9  +  =  cos  cob  <t>  +  sin  ^  sin  ^  cos  A' 

(65) 


Now 


and 


AG 


A*  =  +  AA° 

c  z 


where  Az,  AA  are  both  small. 
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Equation  (65)  becomes 

^2  Z  O 

sin  A0  =  sin  ^  ^  cos  ^  +  cos  ^  ^  ^  sin  ^  sin  AA^  (66) 
and  for  small  quantities 


A9  =  cos  <f>  +  AA°  sin  ^  (67) 


The  crossrange  miss,  MC,  is 

RA0  so 

M^Az 

A  R 

=  r-TK  ^ 

MCa^o 

=  +  R  A  A°  sin  0 
z  ^ 

(68) 

Z 

The  downrange  miss  can  be  obtained  by  solving  for  AX.  The  sine 
law  gives 


sin  (X  -f  AX)  j  H 

■sin(5  +  A0) 

and,  consequently,  keeping  only  first  order  terms  in  AX  and  noting  that 
X  s  ^  gives 

cos  A0  [sin  ^  +  AX  cos  =  sin  cos  AA  (70) 

From  (70)  if  only  first  order  terms  in  AA  and  A0  are  kept,  AX  =  0. 
Therefore,  AX,  the  downrange  miss  attributable  to  Az  and  6A°,  is  of 
second  order  in  these  uncertainties  and  may  thus  be  ignored.  Table  3 
summarizes  the  results  of  the  Phase  II  analysis. 

C.  Numerical  Analysis  of  G8tG  Uncertainties 

The  results  tabulated  in  Tables  2  and  3  could  be  combined  to  give 
analytical  representations  for  target  miss  per  G&G  uncertainty.  The 
results  of  such  a  recombination  would  be  sufficiently  complicated  as  to 
preclude  any  useful  understanding  of  the  relationship  between  GfaG  effects 


( 
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at  the  launch  area  and  target  miss.  Consequently,  the  results  of  the 
previous  sections  were  used  to  generate  curves  giving  the  value  of  target 
miss/GScG  uncertainty  over  wide  limits  of  burnout  angle  and  missile  range. 

The  expressions  in  Tables  2  and  3  are  functions  not  only  of  the 
burnout  angle,  P,  and  range,  ^  but  also  of  h,  i(r,  v,  and  T:  burnout  height, 
downrange  burnout  geocentric  angle,  burnout  velocity,  and  total  time  of 
powered  flight.  Inclusion  of  these  latter  variables  as  independent  variables 
in  the  numerical  analysis  would  increase  the  extent  of  the  required  analysis 
to  a  point  where  any  generality  would  be  lost.  Accordingly,  an  attempt 
was  made  to  specify  these  variables  either  as  constants  or  as  functions  of 
P  and  For  the  other  variables  h,  !(t,  and  T,a  series  of  sample  ICBM 
trajectories  were  examined.  The  time,  T,  was  found  to  weakly  vary  over 
la  rge  variations  in  range  and  burnout  angle  and  consequently  a  constant 
value,  T  =  185  sec,  was  utilized  for  this  analysis.  Sinch  h  and  if  were 
found  to  vary  over  ranges  of  500,  000  -  1,  500,  000  ft  and  4-5  deg,  respec¬ 
tively,  power  series  expansions  were  made  considering  h,  if  to  be  functions 
of  9  and  <t>.  The  series  expansions  were  then  made  to  conform  exactly  to  six 
different  sample  ICBM  trajectories,  three  each  at  5500  nautical  miles  and 
7500  nautical  miles;  each  trajectory  had  a  different  burnout  angle.  These 
trajectories  were  felt  to  be  representative  of  most  operational  ICBM  tra¬ 
jectories  of  interest.  The  resulting  expansions  in  h  and  if  were  used  for 
those  cases  in  which  the  expansions  gave  values  within  the, essentially, 
arbitrary  bounds  of  500,  000  ft  s  h  ^  1, 500,  000  ft,  1*^  s  ijr  s  10°;  if  the  value 
of  h  or  if  exceeded  these  bounds,  the  value  used  in  the  computation  was 
the  relevant  upper  or  lower  bound.  These  more  or  less  artificial  con¬ 
straints  were  found  to  be  necessary  for  those  values  of  B  and  p  widely 
removed  from  the  region  in  which  the  expansions  were  defined. 

For  p  >  TT  (R  >  10,  800  naut  mi)  the  lower  bound  for  h .  500,  000 

feet . was  used  entirely,  since  ballistic  systems  utilized  for  these 

ranges  would  be  under  severe  propulsion  requirements  and,  therefore,  the 

jje 

From  the  "hit  equation,"  (1),  burnout  velocity  can  immediately  be  deter¬ 
mined  as  a  function  of  8  and  p. 


powered  flight  trajectory  might  be  expected  to  be  shaped,  as  nearly  as 
possible,  horizontally.  For  these  ranges  f  was  set  equal  to  10  degrees 
again  because  the  severe  energy  requirements  would  engender  a  powered 
trajectory  extended  to  the  limit  downrange. 

The  curves  so  generated  are  given  below  in  the  following  order: 

1.  Inertially  guided  systems,  500  -  10,  000  naut  mi 

2.  Inertially  guided  systems,  10,  000  -  22,  000  naut  mi 

3.  Radio  guided  systems,  500  -  10,  000  naut  mi 

4.  Radio  guided  systems,  10,  000  -  22,  000  naut  mi 

Curves  of  equal  burnout  velocity  have  also  been  included  in  the 
plots  listed  above  as  an  aid  in  analysis.  For  ranges  in  excess  of  10,  000 
nautical  miles  certain  portions  of  the  P,  R  plot  are  inadmissiable  for 
ballistic  missiles;  these  regions  are  also  indicated  on  the  appropriate 
contour  plots. 


Index  of  Launch  Area  Uncertainty  Contour  Plots 


Solid  Contours  give  Target  Miss/Uncertainty  as  Specified  in  Table  Below- 
Dashed  Contours  give  the  Burnout  Velocity  Requirement  in  Feet/Second 


Guidance 

Type 

Launch  Area  Uncertainty 

Contours  given  in 
Units  of; 

Contour 
Plot  No. 

500  -  10,  500  Nautical 

Miles 

Inertial 

Deflection  of  the  Vertical  (6||) 

Mdr  (ft)/ 6 11  {s^) 

1 

Inertial 

Elevation  above  Ellipsoid  (e) 

Mdr  (ft)/e  (ft) 

2 

Inertial 

Horizontal  Position  (d) 

(ft)/d  (ft) 

3 

Inertial 

Deflection  of  the  Vertical  (6_l) 

Mcr  (ft)/6j,  i^c) 

4 

Inertial 

Horizontal  Position  (c) 

(ft)/c  (ft) 

5 

Inertial 

Azimuth  (A  ) 

'  z 

Mcr  (ft)/A^  (sec) 

6 

10,  000  -  Z2,  000  Nautical  Miles 


Inertial 

Deflection  of  the  Vertical  (6||) 

Mdr  (ft)/6||  (s^) 

7 

Inertial 

Elevation  above  Ellipsoid  (e) 

Mdr  (ft)/e  (ft) 

8 

Ine  rtial 

Horizontal  Position  (d) 

Mdj^  (ft)/d  (ft) 

9 

Inertial 

Deflection  of  the  Vertical  (6j_) 

Mcr  (ft)/6_L  (s^) 

10 

Inertial 

Horizontal  Position  (c) 

(ft)/c  (ft) 

11 

Inertial 

Azimuth  (A  ) 
z 

Mcr  (ft)/A^  (s^) 

12 
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Guidance 

Type 

Launch  Area  Uncertainty 

Contours  given  in 
Units  of: 

Contour 
Plot  No, 

500  -  10,  500  Nautical  Miles 

Radio 

Deflection  of  the  Vertical  (6||) 

Mdr  (ft)/S||  (8^) 

13 

Radio 

Elevation  above  Ellipsoid  (e) 

Mdr  (ft)/e  (ft) 

14 

Radio 

Horizontal  Position  (d) 

(ft)/d  (ft) 

15 

Radio 

Deflection  of  the  Vertical  (6_l) 

Mqr  (ft)/6_L  (sec) 

16 

Radio 

Horizontal  Position  (c) 

Mcr  (ft)/c  (ft) 

17 

Radio 

Azimuth  (A  ) 

'  z 

Mcr  {ft)/A^  (sec) 

18 

10,  000  -  22,  000  Nautical  Miles 

Radio 

Deflection  of  the  Vertical  (6||) 

Mdr  (ft)/ 6 II  (^) 

19 

Radio 

Elevation  above  Ellipsoid  (G) 

Mdr  (ft)/e  (ft) 

20 

Radio 

Horizontal  Position  (d) 

Mjjj^  (ft)/d  (ft) 

21 

Radio 

Deflection  of  the  Vertical  (6_l) 

Mcr  (ft)/6j^  (^ec) 

22 

Radio 

Horizontal  Position  (c) 

Mcr  (ft)/c  (ft) 

23 

Radio 

Azimuth  (A  ) 
z 

Mcr  (ft)/A^  (s^) 

24 

ijinr 
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